3.3.89 \(\int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx\) [289]

3.3.89.1 Optimal result
3.3.89.2 Mathematica [C] (warning: unable to verify)
3.3.89.3 Rubi [A] (verified)
3.3.89.4 Maple [C] (warning: unable to verify)
3.3.89.5 Fricas [C] (verification not implemented)
3.3.89.6 Sympy [F]
3.3.89.7 Maxima [F]
3.3.89.8 Giac [F]
3.3.89.9 Mupad [F(-1)]

3.3.89.1 Optimal result

Integrand size = 25, antiderivative size = 154 \[ \int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx=\frac {2 a^2 \arctan \left (\sqrt {\sin (c+d x)}\right ) \sqrt {e \csc (c+d x)} \sqrt {\sin (c+d x)}}{d}+\frac {2 a^2 \text {arctanh}\left (\sqrt {\sin (c+d x)}\right ) \sqrt {e \csc (c+d x)} \sqrt {\sin (c+d x)}}{d}+\frac {3 a^2 \sqrt {e \csc (c+d x)} \operatorname {EllipticF}\left (\frac {1}{2} \left (c-\frac {\pi }{2}+d x\right ),2\right ) \sqrt {\sin (c+d x)}}{d}+\frac {a^2 \sqrt {e \csc (c+d x)} \tan (c+d x)}{d} \]

output
2*a^2*arctan(sin(d*x+c)^(1/2))*(e*csc(d*x+c))^(1/2)*sin(d*x+c)^(1/2)/d+2*a 
^2*arctanh(sin(d*x+c)^(1/2))*(e*csc(d*x+c))^(1/2)*sin(d*x+c)^(1/2)/d-3*a^2 
*(sin(1/2*c+1/4*Pi+1/2*d*x)^2)^(1/2)/sin(1/2*c+1/4*Pi+1/2*d*x)*EllipticF(c 
os(1/2*c+1/4*Pi+1/2*d*x),2^(1/2))*(e*csc(d*x+c))^(1/2)*sin(d*x+c)^(1/2)/d+ 
a^2*(e*csc(d*x+c))^(1/2)*tan(d*x+c)/d
 
3.3.89.2 Mathematica [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 5 vs. order 4 in optimal.

Time = 12.89 (sec) , antiderivative size = 168, normalized size of antiderivative = 1.09 \[ \int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx=-\frac {2 a^2 \cos ^5\left (\frac {1}{2} (c+d x)\right ) \sqrt {e \csc (c+d x)} \left (-1+2 \arctan \left (\sqrt {\csc (c+d x)}\right ) \sqrt {\cos ^2(c+d x)} \sqrt {\csc (c+d x)}-2 \text {arctanh}\left (\sqrt {\csc (c+d x)}\right ) \sqrt {\cos ^2(c+d x)} \sqrt {\csc (c+d x)}+3 \sqrt {-\cot ^2(c+d x)} \operatorname {Hypergeometric2F1}\left (\frac {1}{4},\frac {1}{2},\frac {5}{4},\csc ^2(c+d x)\right )\right ) \sec (c+d x) \sec ^4\left (\frac {1}{2} \csc ^{-1}(\csc (c+d x))\right ) \sin \left (\frac {1}{2} (c+d x)\right )}{d} \]

input
Integrate[Sqrt[e*Csc[c + d*x]]*(a + a*Sec[c + d*x])^2,x]
 
output
(-2*a^2*Cos[(c + d*x)/2]^5*Sqrt[e*Csc[c + d*x]]*(-1 + 2*ArcTan[Sqrt[Csc[c 
+ d*x]]]*Sqrt[Cos[c + d*x]^2]*Sqrt[Csc[c + d*x]] - 2*ArcTanh[Sqrt[Csc[c + 
d*x]]]*Sqrt[Cos[c + d*x]^2]*Sqrt[Csc[c + d*x]] + 3*Sqrt[-Cot[c + d*x]^2]*H 
ypergeometric2F1[1/4, 1/2, 5/4, Csc[c + d*x]^2])*Sec[c + d*x]*Sec[ArcCsc[C 
sc[c + d*x]]/2]^4*Sin[(c + d*x)/2])/d
 
3.3.89.3 Rubi [A] (verified)

Time = 0.56 (sec) , antiderivative size = 109, normalized size of antiderivative = 0.71, number of steps used = 7, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.280, Rules used = {3042, 4366, 3042, 4360, 3042, 3352, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (a \sec (c+d x)+a)^2 \sqrt {e \csc (c+d x)} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \left (a-a \csc \left (c+d x-\frac {\pi }{2}\right )\right )^2 \sqrt {e \sec \left (c+d x-\frac {\pi }{2}\right )}dx\)

\(\Big \downarrow \) 4366

\(\displaystyle \sqrt {\sin (c+d x)} \sqrt {e \csc (c+d x)} \int \frac {(\sec (c+d x) a+a)^2}{\sqrt {\sin (c+d x)}}dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \sqrt {\sin (c+d x)} \sqrt {e \csc (c+d x)} \int \frac {\left (a-a \csc \left (c+d x-\frac {\pi }{2}\right )\right )^2}{\sqrt {\cos \left (c+d x-\frac {\pi }{2}\right )}}dx\)

\(\Big \downarrow \) 4360

\(\displaystyle \sqrt {\sin (c+d x)} \sqrt {e \csc (c+d x)} \int \frac {(-\cos (c+d x) a-a)^2 \sec ^2(c+d x)}{\sqrt {\sin (c+d x)}}dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \sqrt {\sin (c+d x)} \sqrt {e \csc (c+d x)} \int \frac {\left (a \sin \left (c+d x-\frac {\pi }{2}\right )-a\right )^2}{\sqrt {\cos \left (c+d x-\frac {\pi }{2}\right )} \sin \left (c+d x-\frac {\pi }{2}\right )^2}dx\)

\(\Big \downarrow \) 3352

\(\displaystyle \sqrt {\sin (c+d x)} \sqrt {e \csc (c+d x)} \int \left (\frac {\sec ^2(c+d x) a^2}{\sqrt {\sin (c+d x)}}+\frac {2 \sec (c+d x) a^2}{\sqrt {\sin (c+d x)}}+\frac {a^2}{\sqrt {\sin (c+d x)}}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \sqrt {\sin (c+d x)} \sqrt {e \csc (c+d x)} \left (\frac {2 a^2 \arctan \left (\sqrt {\sin (c+d x)}\right )}{d}+\frac {2 a^2 \text {arctanh}\left (\sqrt {\sin (c+d x)}\right )}{d}+\frac {3 a^2 \operatorname {EllipticF}\left (\frac {1}{2} \left (c+d x-\frac {\pi }{2}\right ),2\right )}{d}+\frac {a^2 \sqrt {\sin (c+d x)} \sec (c+d x)}{d}\right )\)

input
Int[Sqrt[e*Csc[c + d*x]]*(a + a*Sec[c + d*x])^2,x]
 
output
Sqrt[e*Csc[c + d*x]]*((2*a^2*ArcTan[Sqrt[Sin[c + d*x]]])/d + (2*a^2*ArcTan 
h[Sqrt[Sin[c + d*x]]])/d + (3*a^2*EllipticF[(c - Pi/2 + d*x)/2, 2])/d + (a 
^2*Sec[c + d*x]*Sqrt[Sin[c + d*x]])/d)*Sqrt[Sin[c + d*x]]
 

3.3.89.3.1 Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3352
Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((d_.)*sin[(e_.) + (f_.)*(x_)])^(n 
_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Int[ExpandTrig 
[(g*cos[e + f*x])^p, (d*sin[e + f*x])^n*(a + b*sin[e + f*x])^m, x], x] /; F 
reeQ[{a, b, d, e, f, g, n, p}, x] && EqQ[a^2 - b^2, 0] && IGtQ[m, 0]
 

rule 4360
Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_.)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + 
(a_))^(m_.), x_Symbol] :> Int[(g*Cos[e + f*x])^p*((b + a*Sin[e + f*x])^m/Si 
n[e + f*x]^m), x] /; FreeQ[{a, b, e, f, g, p}, x] && IntegerQ[m]
 

rule 4366
Int[(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_.)*((g_.)*sec[(e_.) + (f_.)*( 
x_)])^(p_), x_Symbol] :> Simp[g^IntPart[p]*(g*Sec[e + f*x])^FracPart[p]*Cos 
[e + f*x]^FracPart[p]   Int[(a + b*Csc[e + f*x])^m/Cos[e + f*x]^p, x], x] / 
; FreeQ[{a, b, e, f, g, m, p}, x] &&  !IntegerQ[p]
 
3.3.89.4 Maple [C] (warning: unable to verify)

Result contains complex when optimal does not.

Time = 12.56 (sec) , antiderivative size = 476, normalized size of antiderivative = 3.09

method result size
parts \(\frac {i a^{2} \left (\cos \left (d x +c \right )+1\right ) \sqrt {2}\, \sqrt {e \csc \left (d x +c \right )}\, \sqrt {-i \left (i-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )}\, \sqrt {-i \left (i+\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )}\, \sqrt {-i \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )}\, \operatorname {EllipticF}\left (\sqrt {-i \left (i-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )}, \frac {\sqrt {2}}{2}\right )}{d}+\frac {a^{2} \sqrt {2}\, \sqrt {e \csc \left (d x +c \right )}\, \left (i \sqrt {-i \left (i-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )}\, \sqrt {-i \left (i+\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )}\, \sqrt {-i \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )}\, \operatorname {EllipticF}\left (\sqrt {-i \left (i-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )}, \frac {\sqrt {2}}{2}\right ) \cos \left (d x +c \right )+i \sqrt {-i \left (i-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )}\, \sqrt {-i \left (i+\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )}\, \sqrt {-i \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )}\, \operatorname {EllipticF}\left (\sqrt {-i \left (i-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )}, \frac {\sqrt {2}}{2}\right )+\sqrt {2}\, \tan \left (d x +c \right )\right )}{2 d}+\frac {2 a^{2} \sin \left (d x +c \right ) \left (\operatorname {arctanh}\left (\sqrt {\frac {\sin \left (d x +c \right )}{\left (\cos \left (d x +c \right )+1\right )^{2}}}\, \left (\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )\right )-\arctan \left (\sqrt {\frac {\sin \left (d x +c \right )}{\left (\cos \left (d x +c \right )+1\right )^{2}}}\, \left (\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )\right )\right ) \sqrt {e \csc \left (d x +c \right )}}{d \left (\cos \left (d x +c \right )+1\right ) \sqrt {\frac {\sin \left (d x +c \right )}{\left (\cos \left (d x +c \right )+1\right )^{2}}}}\) \(476\)
default \(\text {Expression too large to display}\) \(1034\)

input
int((a+a*sec(d*x+c))^2*(e*csc(d*x+c))^(1/2),x,method=_RETURNVERBOSE)
 
output
I*a^2/d*(cos(d*x+c)+1)*2^(1/2)*(e*csc(d*x+c))^(1/2)*(-I*(I-cot(d*x+c)+csc( 
d*x+c)))^(1/2)*(-I*(I+cot(d*x+c)-csc(d*x+c)))^(1/2)*(-I*(cot(d*x+c)-csc(d* 
x+c)))^(1/2)*EllipticF((-I*(I-cot(d*x+c)+csc(d*x+c)))^(1/2),1/2*2^(1/2))+1 
/2*a^2/d*2^(1/2)*(e*csc(d*x+c))^(1/2)*(I*(-I*(I-cot(d*x+c)+csc(d*x+c)))^(1 
/2)*(-I*(I+cot(d*x+c)-csc(d*x+c)))^(1/2)*(-I*(cot(d*x+c)-csc(d*x+c)))^(1/2 
)*EllipticF((-I*(I-cot(d*x+c)+csc(d*x+c)))^(1/2),1/2*2^(1/2))*cos(d*x+c)+I 
*(-I*(I-cot(d*x+c)+csc(d*x+c)))^(1/2)*(-I*(I+cot(d*x+c)-csc(d*x+c)))^(1/2) 
*(-I*(cot(d*x+c)-csc(d*x+c)))^(1/2)*EllipticF((-I*(I-cot(d*x+c)+csc(d*x+c) 
))^(1/2),1/2*2^(1/2))+2^(1/2)*tan(d*x+c))+2*a^2/d*sin(d*x+c)*(arctanh((sin 
(d*x+c)/(cos(d*x+c)+1)^2)^(1/2)*(cot(d*x+c)+csc(d*x+c)))-arctan((sin(d*x+c 
)/(cos(d*x+c)+1)^2)^(1/2)*(cot(d*x+c)+csc(d*x+c))))*(e*csc(d*x+c))^(1/2)/( 
cos(d*x+c)+1)/(sin(d*x+c)/(cos(d*x+c)+1)^2)^(1/2)
 
3.3.89.5 Fricas [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.20 (sec) , antiderivative size = 664, normalized size of antiderivative = 4.31 \[ \int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx=\left [-\frac {2 \, a^{2} \sqrt {-e} \arctan \left (-\frac {{\left (\cos \left (d x + c\right )^{2} - 6 \, \sin \left (d x + c\right ) - 2\right )} \sqrt {-e} \sqrt {\frac {e}{\sin \left (d x + c\right )}}}{4 \, {\left (e \sin \left (d x + c\right ) + e\right )}}\right ) \cos \left (d x + c\right ) - a^{2} \sqrt {-e} \cos \left (d x + c\right ) \log \left (\frac {e \cos \left (d x + c\right )^{4} - 72 \, e \cos \left (d x + c\right )^{2} + 8 \, {\left (\cos \left (d x + c\right )^{4} - 9 \, \cos \left (d x + c\right )^{2} + {\left (7 \, \cos \left (d x + c\right )^{2} - 8\right )} \sin \left (d x + c\right ) + 8\right )} \sqrt {-e} \sqrt {\frac {e}{\sin \left (d x + c\right )}} + 28 \, {\left (e \cos \left (d x + c\right )^{2} - 2 \, e\right )} \sin \left (d x + c\right ) + 72 \, e}{\cos \left (d x + c\right )^{4} - 8 \, \cos \left (d x + c\right )^{2} - 4 \, {\left (\cos \left (d x + c\right )^{2} - 2\right )} \sin \left (d x + c\right ) + 8}\right ) + 6 i \, a^{2} \sqrt {2 i \, e} \cos \left (d x + c\right ) {\rm weierstrassPInverse}\left (4, 0, \cos \left (d x + c\right ) + i \, \sin \left (d x + c\right )\right ) - 6 i \, a^{2} \sqrt {-2 i \, e} \cos \left (d x + c\right ) {\rm weierstrassPInverse}\left (4, 0, \cos \left (d x + c\right ) - i \, \sin \left (d x + c\right )\right ) - 4 \, a^{2} \sqrt {\frac {e}{\sin \left (d x + c\right )}} \sin \left (d x + c\right )}{4 \, d \cos \left (d x + c\right )}, -\frac {2 \, a^{2} \sqrt {e} \arctan \left (\frac {{\left (\cos \left (d x + c\right )^{2} + 6 \, \sin \left (d x + c\right ) - 2\right )} \sqrt {e} \sqrt {\frac {e}{\sin \left (d x + c\right )}}}{4 \, {\left (e \sin \left (d x + c\right ) - e\right )}}\right ) \cos \left (d x + c\right ) - a^{2} \sqrt {e} \cos \left (d x + c\right ) \log \left (\frac {e \cos \left (d x + c\right )^{4} - 72 \, e \cos \left (d x + c\right )^{2} + 8 \, {\left (\cos \left (d x + c\right )^{4} - 9 \, \cos \left (d x + c\right )^{2} - {\left (7 \, \cos \left (d x + c\right )^{2} - 8\right )} \sin \left (d x + c\right ) + 8\right )} \sqrt {e} \sqrt {\frac {e}{\sin \left (d x + c\right )}} - 28 \, {\left (e \cos \left (d x + c\right )^{2} - 2 \, e\right )} \sin \left (d x + c\right ) + 72 \, e}{\cos \left (d x + c\right )^{4} - 8 \, \cos \left (d x + c\right )^{2} + 4 \, {\left (\cos \left (d x + c\right )^{2} - 2\right )} \sin \left (d x + c\right ) + 8}\right ) + 6 i \, a^{2} \sqrt {2 i \, e} \cos \left (d x + c\right ) {\rm weierstrassPInverse}\left (4, 0, \cos \left (d x + c\right ) + i \, \sin \left (d x + c\right )\right ) - 6 i \, a^{2} \sqrt {-2 i \, e} \cos \left (d x + c\right ) {\rm weierstrassPInverse}\left (4, 0, \cos \left (d x + c\right ) - i \, \sin \left (d x + c\right )\right ) - 4 \, a^{2} \sqrt {\frac {e}{\sin \left (d x + c\right )}} \sin \left (d x + c\right )}{4 \, d \cos \left (d x + c\right )}\right ] \]

input
integrate((a+a*sec(d*x+c))^2*(e*csc(d*x+c))^(1/2),x, algorithm="fricas")
 
output
[-1/4*(2*a^2*sqrt(-e)*arctan(-1/4*(cos(d*x + c)^2 - 6*sin(d*x + c) - 2)*sq 
rt(-e)*sqrt(e/sin(d*x + c))/(e*sin(d*x + c) + e))*cos(d*x + c) - a^2*sqrt( 
-e)*cos(d*x + c)*log((e*cos(d*x + c)^4 - 72*e*cos(d*x + c)^2 + 8*(cos(d*x 
+ c)^4 - 9*cos(d*x + c)^2 + (7*cos(d*x + c)^2 - 8)*sin(d*x + c) + 8)*sqrt( 
-e)*sqrt(e/sin(d*x + c)) + 28*(e*cos(d*x + c)^2 - 2*e)*sin(d*x + c) + 72*e 
)/(cos(d*x + c)^4 - 8*cos(d*x + c)^2 - 4*(cos(d*x + c)^2 - 2)*sin(d*x + c) 
 + 8)) + 6*I*a^2*sqrt(2*I*e)*cos(d*x + c)*weierstrassPInverse(4, 0, cos(d* 
x + c) + I*sin(d*x + c)) - 6*I*a^2*sqrt(-2*I*e)*cos(d*x + c)*weierstrassPI 
nverse(4, 0, cos(d*x + c) - I*sin(d*x + c)) - 4*a^2*sqrt(e/sin(d*x + c))*s 
in(d*x + c))/(d*cos(d*x + c)), -1/4*(2*a^2*sqrt(e)*arctan(1/4*(cos(d*x + c 
)^2 + 6*sin(d*x + c) - 2)*sqrt(e)*sqrt(e/sin(d*x + c))/(e*sin(d*x + c) - e 
))*cos(d*x + c) - a^2*sqrt(e)*cos(d*x + c)*log((e*cos(d*x + c)^4 - 72*e*co 
s(d*x + c)^2 + 8*(cos(d*x + c)^4 - 9*cos(d*x + c)^2 - (7*cos(d*x + c)^2 - 
8)*sin(d*x + c) + 8)*sqrt(e)*sqrt(e/sin(d*x + c)) - 28*(e*cos(d*x + c)^2 - 
 2*e)*sin(d*x + c) + 72*e)/(cos(d*x + c)^4 - 8*cos(d*x + c)^2 + 4*(cos(d*x 
 + c)^2 - 2)*sin(d*x + c) + 8)) + 6*I*a^2*sqrt(2*I*e)*cos(d*x + c)*weierst 
rassPInverse(4, 0, cos(d*x + c) + I*sin(d*x + c)) - 6*I*a^2*sqrt(-2*I*e)*c 
os(d*x + c)*weierstrassPInverse(4, 0, cos(d*x + c) - I*sin(d*x + c)) - 4*a 
^2*sqrt(e/sin(d*x + c))*sin(d*x + c))/(d*cos(d*x + c))]
 
3.3.89.6 Sympy [F]

\[ \int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx=a^{2} \left (\int \sqrt {e \csc {\left (c + d x \right )}}\, dx + \int 2 \sqrt {e \csc {\left (c + d x \right )}} \sec {\left (c + d x \right )}\, dx + \int \sqrt {e \csc {\left (c + d x \right )}} \sec ^{2}{\left (c + d x \right )}\, dx\right ) \]

input
integrate((a+a*sec(d*x+c))**2*(e*csc(d*x+c))**(1/2),x)
 
output
a**2*(Integral(sqrt(e*csc(c + d*x)), x) + Integral(2*sqrt(e*csc(c + d*x))* 
sec(c + d*x), x) + Integral(sqrt(e*csc(c + d*x))*sec(c + d*x)**2, x))
 
3.3.89.7 Maxima [F]

\[ \int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx=\int { \sqrt {e \csc \left (d x + c\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{2} \,d x } \]

input
integrate((a+a*sec(d*x+c))^2*(e*csc(d*x+c))^(1/2),x, algorithm="maxima")
 
output
integrate(sqrt(e*csc(d*x + c))*(a*sec(d*x + c) + a)^2, x)
 
3.3.89.8 Giac [F]

\[ \int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx=\int { \sqrt {e \csc \left (d x + c\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{2} \,d x } \]

input
integrate((a+a*sec(d*x+c))^2*(e*csc(d*x+c))^(1/2),x, algorithm="giac")
 
output
integrate(sqrt(e*csc(d*x + c))*(a*sec(d*x + c) + a)^2, x)
 
3.3.89.9 Mupad [F(-1)]

Timed out. \[ \int \sqrt {e \csc (c+d x)} (a+a \sec (c+d x))^2 \, dx=\int {\left (a+\frac {a}{\cos \left (c+d\,x\right )}\right )}^2\,\sqrt {\frac {e}{\sin \left (c+d\,x\right )}} \,d x \]

input
int((a + a/cos(c + d*x))^2*(e/sin(c + d*x))^(1/2),x)
 
output
int((a + a/cos(c + d*x))^2*(e/sin(c + d*x))^(1/2), x)